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Abstract 

Leading corrections to Planck's formula and photon thermodynamics aris- 
ing from the pair-mediated photon-photon interaction are calculated. This 
interaction is attractive and causes an increase in occupation number for all 
modes. Possible consequences, including the role of the cosmic photon gas in 
structure formation, are considered. [lhep-th/9308089|| 
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I. INTRODUCTION 



Electromagnetic radiation in thermal equilibrium is a fundamental physical system which 
has played a crucial, if serendipitous, role in opening up new frontiers in physics; few systems 
are as ubiquitous or have played as significant a role in the progress of modern physics 
and astronomy. Indeed it was in an attempt to describe the properties of cavity radiation 
that Planck stumbled upon the notion of quantization in 1900, and, 65 years later, it was 
while testing a microwave antenna designed for communications research that Penzias and 
Wilson discovered the cosmic microwave background radiation, the single most important 
clue to the past history of the universe. The recent discovery by Smoot and coworkers [1] of 
large-angle temperature fluctuations in the background radiation, which, by contrast, had 
long been anticipated, promises to be another milestone in the quest for the origins of the 
universe. The distribution law discovered by Planck accurately describes the equilibrium 
properties of an assembly of photons over a vast range of temperatures and scales, from 
terrestrial cavity radiation (from which it was deduced) to hot stellar atmospheres, and, 
of course, including the 2.73K cosmic background radiation. Nevertheless, there are small, 
high-temperature corrections to Planck's law and photon thermodynamics which, in view of 
the pervasive presence of radiation systems in general and the singularly important role of 
the background radiation for cosmic evolution in particular, may play a significant role in an 
as yet undiscovered phenomenon. For this reason, as well as their intrinsic interest as high 
temperature modifications to a fundamental law of nature, these corrections are worthy of 
serious attention. This communication is devoted to the calculation of these corrections and 
the exploration of some of their consequences. 

The Planck law basically describes a noninteracting gas of massless, spin-one bosons (not 
subject to number conservation), driven to equilibrium by interaction with an external sys- 
tem, e.g., the atoms in the walls of the cavity in the case of cavity radiation, and the charged 
particles present in the cosmic fluid in the case of the background radiation. The virtual 
nonexistence of an interaction between photons under ordinary conditions, a hallmark of 
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Maxwell electrodynamics and the working principle of virtually all existing telecommunica- 
tion systems, is what sets the photon gas apart from other ideal gases. To highlight this 
difference, let us consider the fact that the mean free path for hydrogen under standard 
conditions is about 10~^ m, allowing hydrogen to behave like a continuous medium on scales 
larger than a micron. By contrast, the same quantity for a photon gas is a staggering 10^^ 
m [cf. Eq. (19) below], implying that under ordinary conditions photons behave like a 
collection of independent, free-streaming particles on any physically meaningful scale. 

The absence of interaction between photons and the linearity of electrodynamics are of 
course violated by the process of pair creation, a fact that was realized by Halpern [2] as 
early as in 1933. Not long before this, Delbriick [3] had suggested that quantum effects 
would cause photons to be scattered by an external electric field. Delbriick's phenomenon, 
which was observed in the scattering of 1.33 Mev gamma rays by the Coulomb field of the 
lead nucleus some years later [4] , is essentially the interaction of a real photon with a virtual 
one [5]. By contrast, the interaction of two real photons predicted by Halpern has never 
been detected because of its extreme weakness. The center-of-mass photon-photon cross 
section is equal to .03lQ;^rg(a;/m)^ for uj/m <^ 1 and nearly equal to it almost up to the pair 
production threshold a; = m [6]; here a is the fine-structure constant, re the classical electron 
radius, uj the photon energy, and m the electron mass (natural units with h = c = k = 1 will 
be used unless otherwise stated). For 633 nm light, this cross section is about 4 x 10~^^ cm^, 
and it rises steeply to about 3 x 10~^° cm^ near the production threshold. These minute cross 
sections are the reason why phenomena involving the interaction of real photons have not 
been much discussed in connection with nonlinear electrodynamic effects, attention having 
been focused instead on the magnetic counterpart of Delbriick scattering and other effects 
[7]. It is reasonable to expect, however, that these tiny magnitudes will eventually succumb 
to high-precision experiments, or to new observations. Indeed a scheme has already been 
proposed for detecting the photon-photon interaction by colliding laser beam interferometry, 
using a modification of Michclson interferometers developed for gravity wave detection [8]. 

The present discussion deals with photon interactions in equilibrium distributions, the 
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opposite extreme to highly-ordered, far-from-equihbrium distributions characteristic of laser 
light. While at low (i.e., -C m) temperatures a photon gas is basically pure, at sufficiently 
high temperatures a significant contamination of electron-positron (and possibly other) pairs 
may be present, making it inappropriate to treat the gas as a single-component system. To 
ensure that the photon gas is effectively free of contamination, we shall restrict our treatment 
to temperatures for which the photon-photon collision rate dominates the collision rate of the 
photons with the pair-produced electrons or positrons so that the latter play no significant 
role in the thermodynamics of the photon gas. We shall see below that the condition 
T < 1.5 X 10® K, where T is the temperature of the photon gas, is sufficient for this purpose. 
Now for this range of temperatures, the parameter T/m is no larger than 2.5 x 10~^ and 
can be treated as small. Since, for dimensional reasons, the dependence of the corrections 
we are seeking on T must be through the combination T/m, we can achieve a considerable 
simplification in the calculations by limiting them to the leading order in this parameter, with 
essentially no loss in accuracy. Our objective, therefore, is to derive the thermodynamics 
of a photon gas (i.e., the charge-zero, equilibrium state of QED) up to and including the 
leading correction in the parameters a and T/m. 

The leading contribution to photon-photon scattering is of second order in a and orig- 
inates in the box diagram in which four external photon lines are attached to a closed 
electron loop [9]. In leading order, the matrix element for this diagram is directly related 
to the interaction energy of a pair of photons, so that the corrections in question will be 
of relative order Q;^(T/m)^, the power of 4 resulting from the four electron lines of the box 
diagram in the limit T/m 0. Obviously, these corrections will only be manifested un- 
der extreme conditions, requiring very high temperatures or very large spatial dimensions. 
Such conditions may be achievable in controlled experiments, or may have prevailed during 
the evolution of the universe subsequent to the annihilation of positrons. For example, the 
elastic properties of the cosmic photon gas may have played a role in the early stages of 
density perturbation growth and structure formation . As will be shown below, the cosmic 
background radiation has in the past been under conditions such that the photon-photon 
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collision rate exceeded the expansion rate. This fact implies the possibility of independently 
propagating sound waves or supporting growing density perturbations for the photon gas. 
Accordingly, the results derived below will be explored in some detail with a view to their 
cosmological implications, although there do not appear to be any observable consequences 
at this time. 



The thermodynamics of a photon gas can be conveniently derived in the canonical en- 
semble where the density matrix and the partition function are given by 



Here V is the volume and /3 — T~^. To the accuracy sought here, the Hamiltonian may 
be written ets H — Hq + Hint, where Hq describes noninter acting photons and Hint is the 
effective interaction Hamiltonian resulting from the box diagram. Recall that our objective 
is to calculate p in the limit T/m ^ 0. Since the amplitude resulting from the box diagram is 
convergent, this limit can be implemented by taking the static limit, m — > oo, of Hint in Eqs. 
(1) and (2), a step that considerably simplifies the calculation. One simply calculates the 
amplitude for the box diagram in the static limit, then converts this into Hint, in much the 
same way that the Born amplitude is converted into an interaction potential when dealing 
with potential scattering. 

A convenient alternative to the above procedure is to take advantage of a well known 
result, known as the Heisenberg-Euler Lagrangian [10], that recasts the static hmit of the 
box diagram amplitude as an effective interaction Lagrangian expressed in terms of the field 
variables. The density corresponding to this Lagrangian is 



II. PHOTON-PHOTON INTERACTION ENERGY 



pi(3,V) ^ Z-\(3,V)eM-mV)], 



(1) 



Z{P,V)^tr{eM-PH{V)]}. 



(2) 




(3) 



5 



A simple calculation shows that with C^nt treated to first order, the interaction Hamiltonian 
density is given by Hint = —Ant, i-e., ii is equal to |(n^ + B^) — £int(n, B). Note that the 
canonical momentum density H is now equal to E + (?Ant(E, B)/(?E, the second term being 
the effect of the interaction [11]. 

The next step is to calculate the diagonal matrix elements of the interaction Hamiltonian 
Hint in the photon number representation (these being the only ones needed in the following 
calculations). It is convenient to carry out this calculation in the Coulomb gauge using the 
complex field G = 11 + iB. The plane- wave expansion of G in finite- volume normalization 
appears as 

G(i,x) = V ^^^fA(k) 

X [«k,A(^) exp(ik • x) - dl^xi't) exp(-ik • x) , 

where k = |k| and k = k/k. Here the unit vectors 6^(6;), A = 1, 2, represent the two polar- 
ization directions normal to k, and the complex vector fx(k) is given by e;^(k) -|- ikxeA(k). 
Furthermore, ak,A ^-nd ay^ ^ are standard destruction and creation operators subject to the 
commutation relations 



It is useful to note here that the product fA(k) • fA'(k') vanishes for k = k'. 

We now proceed to calculate the diagonal matrix elements of (G^)^ and G^(G''^)^. Let 
\n) represent the state with n(k. A) photons of momentum k and polarization eA(k). Then, 
following a lengthy but essentially straightforward calculation, we find 

(n|(G)2(Gt)» = l^Y.i,2^i^2\ii ■ f2|Vn2, 

where the sum runs over all possible values of (ki,Ai) and (k2,A2), and the abbreviations 
fj = fAi(ki) and rij = n{\ii, Aj) have been introduced. The corresponding matrix element for 
(G^)^ is obtained from this formula by replacing |fi • f2|^ with (fi • f2)^. Note that in these 
sums all terms with ki = k2 vanish as a result of the vanishing property of fi • f2 mentioned 
above. 
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The photon-photon interaction energy can now be written down using the above infor- 
mation; the result is 

{n\H,nt\n) = -^^^^^_^r^l,r^2A;lfc2(47^2 + 7X2), (4) 

where TZ and I are the real and imaginary parts of the quantity fi-f2. Note that (n|i7int|n) 
represents a sum of negative two-body interaction terms; i.e., photons attract photons. Fur- 
thermore, both TZ and I vanish for ki = k2 as a result of the vanishing of fi • f2 mentioned 
above; i.e., parallel photons don't interact. The latter property can be traced to the fact that 
the center-of-mass energy of a pair of parallel photons vanishes, implying the same for the 
corresponding cross section. 



III. PHOTON THERMODYNAMICS 

To first order in the interaction, the density matrix of Eq. (1) can be represented as 

p = /' deexp(-e/3^o)(l - /^iYint) exp[(^ - l)(3Ho], (5) 
Jo 

from which the partition function can be calculated using the result given in Eq. (4); 

^ = E{„}(1 - P{ri\Hint\n)) exp(-/3(n|^o|n)), (6) 

where {n\Ho\n) is the energy of an assembly of noninteracting photons, given by the familiar 
expression X^n(k, X)k. The symbol {n} in Eq. (6) indicates an unrestricted summation (since 
photon number is not fixed) over n(k. A) for each possible value of (k. A). 

Thermodynamic quantities can now be calculated as ensemble averages using Eqs. (5) 
and (6). For example, the mean energy U{P,V), defined by tr{Hp), can be calculated by 
applying the formula —{d/d0)\nZ. To find U, we first substitute (4) in (6) and perform 
the sum over {n} to get 



^0, (7) 
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where g{l, 2) is the quantity multiplying the photon numbers in Eq. (4), rio(/c) is the mean 
occupation number for noninteracting photons as given by the Planck formula [exp(/3/c) — 
and Zq is the associated partition function n[l ~ exp(— The product in the 
last expression runs over all possible values of k. 

To complete the calculation of U, we find it convenient to replace ^'(1, 2) in (7) by its aver- 
age, g{ki, k2), over the two polarization directions as well as the orientations of the momenta 
(a permissible operation since uq only depends on k); the result is —22a'^kik2/lS5Vm'^. Us- 
ing this result in (7) and carrying out the remaining sums (in the large-volume, continuum 
limit), we find 



1 + 



Zo, (8) 



135m4 

where uo is the mean energy density of noninteracting photons, equal to 7r^/15/3^. 
The mean energy density u = U/V can now be calculated from (8); the result is 



u 



1 + ^02^^^/^^) 



Uo, (9) 



where T^, — m is equal to 5.9 GK. As expected, the post-Planckian correction in Eq. (9) is 
very small at terrestrially occurring temperatures, except possibly those achievable in exper- 
iments on thermonuclear research. The sign of the correction is also noteworthy. Although 
the interaction is attractive, at a given temperature the interacting photon gas has a higher 
energy density than the noninteracting one. To elucidate this result, let us consider the 
relation 

u^Uo + V-\l - Td/dT)H;nt, (10) 

which is valid to first order in the correction. Here Hint is the mean value of the interaction 
Hamiltonian tr{Hintp) in leading order. A straightforward calculation now yields 

2025~^^ ^ 

Equation (10) shows that if, as is the case here, Hint depends on the temperature more 
strongly than linear, then u will exceed uo in case the interaction is attractive. We must 
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therefore conclude that the negative photon-photon interaction energy is in the present 
instance more than compensated by an increase in the average number of photons in each 
mode, resulting in a net increase in the internal energy. 

To verify the increase in photon occupation number just predicted, we proceed to calcu- 
late the ensemble average of the photon number operator, tr[n(k, \)p]. To the leading order 
we find, after some algebra. 



n{k) = |l + [1 + Mk)] (T/Te)4 no(fc). 



(12) 



2025m 

This formula gives the leading correction to Planck's radiation law arising from vacuum 
polarization. As predicted above, the correction is an increase in occupation number for all 
modes. Note that the increase is more pronounced for higher energy modes, reflecting the 
increase with energy of the photon-photon scattering cross section. 

Other quantities of interest are calculated in a similar manner. For example, the pressure 
for the photon gas is found to be 



P 



1 + 



675 



Po, 



(13) 



where p(po) is the pressure for the interacting (noninteracting) case. Likewise, the speed of 
sound, given by dp/du in the present case, is found with the help of Eqs. (9) and (13). The 
result is 



2025 



(T/Te)^ 



c 

73' 



(14) 



where c is the speed of light in vacuum. The speed of light in the photon gas, on the other 
hand, is found by first calculating the interaction energy cui^t of an "external" photon of 
momentum k with the photon gas using Eq. (4). The result is 

447r2Q;2 



2025 



{T/nTk. 



(15) 



The speed of light in the photon gas, v^, is then found from the formula (1 -|- dujmt/ dk)c; 



1 - 



2025 



c. 



(16) 
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This result implies that a photon gas in equilibrium acts as a linear, isotropic refractive 
medium for the propagation of electromagnetic waves, its index of refraction being equal to 

As already stated, these post-Planckian corrections are very small, even at high tem- 
peratures. For example, at a temperature of lO^^Tg, which is of the order of the highest 
temperatures presently achievable in thermonuclear research, the fractional correction to 
Planck's formula is of the order of 10~^ (one must also remember the contaminations arising 
from the presence of a plasma and the complications caused by the short-lived nature of the 
temperature peak in an actual experiment). Furthermore, as discussed earlier, the above 
results should be supplemented by the condition that requires a photon to collide much more 
frequently with another photon than with the equilibrium population of pair-produced elec- 
trons or positrons. To find the inequality that enforces this condition, we must first calculate 
the photon-photon collision frequency. 

Using the standard result for the number of coUisions in a gas, we find for the coUision 
frequency of a given photon, 

1^77(^1) = V^"'E2^o(A;2)|ki - k2|a(ki,k2), (17) 

where (j(ki, k2) is the scattering cross section for a pair of photons of momenta ki and k2. 

Transforming this expression to the center-of-mass system and using the value of a given 
earlier, we find 

'-^^^ = ^5562^'^^ ^^'^^^ ^^^^ 
for the collision frequency of a photon of momentum k with the other photons. 

The mean collision frequency can now be found by averaging i/^^{k). The result, stated 
as a mean free path, is given by 

X^Jm] = 1.8 X 10-^ { ^^^^ \ . (19) 

This equation gives the mean free path (expressed in meters) in terms of the temperature 
(expressed in Kelvins). Now the condition that ensures a pure photon gas requires that 
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A-y-y <^ A-ye, wheiB A-ye Is thc (partial) mean free path for the coUisions of the photons with 
pair-produced electrons or positrons. 

A calculation involving the equilibrium concentration of electron-positron pairs and 
the Thomson cross section shows that for T/Te <S 1, the ratio A^^/A^e is given by 
A{TjTy^/'^ ex.p{-TjT), where A ~ 5 x 10^. For T = 1.5 x 10^ K, which is the upper 
limit to the range of temperatures considered in this paper, this ratio is less than 4%, so 
that the purity condition is well satisfied for the designated range of temperatures. Thus 
for T < 1.5 X 10^ K, we have in Eqs. (8)-(16) the desired post-Planckian corrections to the 
thermodynamics of an assembly of photons in equilibrium. 

IV. IMPLICATIONS FOR THE COSMIC PHOTON GAS 

At this point the question of the meaning of V in the above derivations must be considered 
more carefully since the naive model of an enclosure with impenetrable walls is clearly not 
tenable at the range of temperatures where post-Planckian corrections become significant. 
Indeed it is necessary to examine the role of whatever mechanism serves to confine the 
photon gas to ensure that it does not appreciably affect the thermodynamics, just as was 
done with pair-produced electrons and positrons. As an example, consider the cosmic photon 
gas, where one might impose the requirement that the scale of inhomogeneities in the cosmic 
gravitational field, which is of the order of the cosmic scale parameter R, be much larger 
than the photon mean free path. This requirement basically translates into the inequality 
A^^ <S y^/^, where F ~ i?^, a condition that can be interpreted to mean that photon- photon 
collisions must be far more frequent than photon-graviton collisions [12]. Using Eq. (19), 
we can rewrite the last condition in terms of the temperature and arrive at the requirements 

V^/^[m\ > 4.0 X 10' 

Observe that these restrictions on temperature and volume also guarantee that on scales 
larger than A^-y the photon gas is in the hydrodynamic regime (in contrast to the independent- 
particle regime prevalent under ordinary conditions) , and essentially behaves like an elastic 
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T[K] 
1.5 X 108 



— 1> 



T < 1.5 X 10^ K. 



(20) 



continuum. For example, such a gas can support or propagate density perturbations of 
spatial dimensions d, provided V^^^ ~:$> d ^ X^^. It is clear, however, that the necessary 
conditions for such phenomena are extreme, suggesting an examination of the early stages 
of the universe. 

According to the standard model [13], about one hour after the big bang the contents of 
the universe had cooled to a temperature of about 1.5 x 10^ K, and the spatio-temporal scales 
of the universe, which for the present discussion may be taken to be the Hubble distance and 
time, were of the order of 10^^ m and 10^ s respectively. The requirements expressed in Eq. 
(20) were thus satisfied. Disregarding for a moment the rest of the contents of the universe 
at this time, we can see that the photon gas behaves hydrodynamically on scales of the 
order of d, provided 10^ m <^ d <S 10^^ m. For example, density perturbations characterized 
by wavelengths and frequencies of the order of 10^° m and 10~^ s~^, respectively, satisfy 
these conditions and could have been supported by the cosmic photon gas. Needless to say, 
the free electrons present in the cosmic fluid at this time (due to the baryon asymmetry of 
the universe) cannot be disregarded since they play a major role by virtue of their frequent 
collisions with the photons. A simple calculation shows that at T ~ 1.5 x 10^ K, the 
ratio u^e/^-y-y is about 10^, and it rapidly increases according to {T/Tg)~^ as the universe 
cools. Indeed by the time the free electrons disappear due to recombination, the photon 
gas is essentially noninteracting, as can be seen by considering Eq. (19). Although at the 
time of T ~ 10^ K the cosmic photon gas was capable of independently supporting density 
perturbations and possibly playing a distinct role in structure formation, its dynamics was 
largely dominated by the free electrons , thereby obviating any interesting effects associated 
with the photons. 

It is amusing in this connection to recall that the phenomenon of sonoluminescence, or 
the creation of light from sound, was observed some sixty years ago. The results of this 
paper suggest that the opposite phenomenon, which may be named the luminosonic effect, 
may also have occurred, albeit early in the life of the universe, and with x-rays as "light" 
and infrasonic waves as "sound." 
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